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ABSTRACT 

This paper considers weighted sum rate maximization of multiuser 
multiple-input single-output interference channel (MISO-IFC) un- 
der outage constraints. The outage-constrained weighted sum rate 
maximization problem is a nonconvex optimization problem and is 
difficult to solve. While it is possible to optimally deal with this 
problem in an exhaustive search manner by finding all the Pareto- 
optimal rate tuples in the (discretized) outage-constrained achiev- 
able rate region, this approach, however, suffers from a prohibitive 
computational complexity and is feasible only when the number of 
transmitter-receive pairs is small. In this paper, we propose a convex 
optimization based approximation method for efficiently handling 
the outage-constrained weighted sum rate maximization problem. 
The proposed approximation method consists of solving a sequence 
of convex optimization problems, and thus can be efficiently imple- 
mented by interior-point methods. Simulation results show that the 
proposed method can yield near-optimal solutions. 

Index Terms — Multiuser interference channel, weighted sum 
rate maximization, outage probability, convex optimization 

1. INTRODUCTION 

Recently, interference management for improving spectral efficiency 
of wireless multiuser systems has been a research topic drawing 
significant attention (TJ. This paper considers the iv-user multiple- 
input single-output interference channel (MISO-IFC) where K 
multi-antenna transmitters simultaneously communicate with K re- 
spective single-antenna receivers over a common frequency band. 
This MISO-IFC arises, for example, in multicell wireless systems 
where each of the base stations is equipped with multiple antennas 
and each mobile station has only one antenna. Under the assumption 
that the transmitters have the perfect channel state information, and 
that the receivers employ single-user detection, it has been shown 
that transmit beamforming is an optimal transmission scheme to 
attain the Pareto boundary of the achievable rate region of MISO- 
IFC (2). The structure of the Pareto-optimal beamforming schemes 
has also been studied in J3] |4j. A game-theoretic approach for 
MISO-IFC has been presented in 0. 

This paper assumes that the channel coefficients are block-faded, 
and that the transmitters know only the statistical distribution of the 
channels. Specifically, each channel is assumed to be circularly 
symmetric complex Gaussian distributed, with a covariance matrix 
known to the transmitters. Under limited delay constraints and due 
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to channel fading, the receivers' performance may suffer from out- 
age. Assuming that the transmitters employ transmit beamforming, 
the achievable rate region of MISO-IFC under outage constraints 
on receivers' performance has been investigated in |6|. While this 
outage-constrained achievable rate region is not known analytically 
so far, it has been shown that this region can be found numerically 
using an exhaustive search method {§]. This method, unfortunately, 
has a complexity that increases exponentially with K{K — 1), and 
therefore is not feasible in practice. 

In this paper, we investigate efficient approaches to achieving 
Pareto-optimal beamforming solutions that maximize the achievable 
weighted sum rate. To this end, we study the design formulation that 
maximizes the weighted sum rate subject to outage constraints and 
individual power constraints. Due to the nonconvextity of the out- 
age constraints, solving the weighted sum rate maximization prob- 
lem is a challenging task. To efficiently deal with this problem, we 
propose a sequential convex approximation method. The proposed 
approximation method is conservative in the sense that the obtained 
approximate beamforming solutions are guaranteed to be feasible 
and satisfy the outage constraints of the original problem. Since the 
proposed method only involves solving convex optimization prob- 
lems, it can be efficiently implemented by interior-point methods in 
a polynomial-time complexity (TJ. The presented simulation results 
show that the proposed approximation method can provide near- 
optimal performance and outperform the existing maximum-ratio 
and zero-forcing transmission strategies. 

2. SIGNAL MODEL AND PROBLEM STATEMENT 

We consider the A'-user MISO interference channel where each of 
the transmitters has Nt antennas and all the receivers are equipped 
with a single antenna. All the transmitters employ transmit beam- 
forming to transmit information signals to their respective receivers. 
Let Si (t) denote the information signal sent from transmitter i, and 
let Wi £ C 4 be the associated beamforming vector. The received 
signal at receiver i is given by 

K 

n(t) = h.fiWiSi(t) + 2^ hkiWkSk(t) + rn(t), (1) 

where hjy £ C^' denotes the channel vector from transmitter k 
to receiver i, and rii(t) is the additive noise of receiver i. The 
noise (t) is assumed to be complex Gaussian distributed with zero 
mean and variance erf > 0, i.e., m(t) ~ CJV(0, of). Assuming 
that Si(t) ~ CA/ r (0, 1) and that the receivers decode the informa- 
tion message using single-user detection (which treats the cross- 



link interference as noise), the achievable rate of the ith transmitter- 
receiver pair is given by 
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In this paper, we assume that the channel coefficients hki are 
block-faded, and that the transmitters can only acquire the statistical 
distribution of the channels. In particular, the elements of hki are 
assumed to be circularly symmetric complex Gaussian distributed 
with covariance matrix equal to Qki b (positive semidefinite), 
i.e., h ki ~ CJV(0, Qki), for all k,i = 1,...,K. Let ft > be 
the target transmission rate of receiver i. Due to channel fading, the 
receivers' performance may suffer from outage; that is, it would have 
a nonzero probability such that n [{hki}^!, {wk}k=i) < Ri- The 
e^-outage achievable rate region is defined as follows: 

Definition 1 [6] Let Pi > denote the power constraint of trans- 
mitter i, and let Ci G (0, 1] denote the maximum tolerable out- 
age probability of receiver i, for i = 1, . . . , K. The rate tuple 
(Ri, ■ ■ ■ , Rk) is said to be achievable if 



Pr jri({h fci }f =1 ,{u>fc}f =1 j <Ri}<€i, i=l,..., 
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for some (wx, . . . , wk) £ Wi X • • • X Wk where Wi = {w G 
C^'j ||iu|| 2 < Pi}- The ei-outage achievable rate region is given by 
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Given the outage specifications ei, . . . ,ejc, it is desirable to 
optimize the beamforming vectors {wk}k=i sucn that the system 
can operate on the so-called Pareto boundary of the achievable rate 
region 1Z (6), with system utilities such as the (weighted) sum of 
Ri, . . . , Rk being maximized. To this end, we consider the follow- 
ing weighted sum rate maximization problem 
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where on > is the priority weight for the ith transmitter-receiver 
pair. Solving problem l[2} is challenging because the outage con- 
straints in d2bb are difficult to handle. One possible approach to solv- 
ing problem (O is to first obtain a set of Pareto-optimal rate tuples 
(Ri, . . . , Rk) by discretizing 1Z using an exhaustive search method 
reported in (6), followed by picking the one that corresponds to the 
largest value of XL = i ctiRi- The complexity of this approach, how- 
ever, increases exponentially with K(K — 1 In the next section, 
based on convex approximation techniques, we present a suboptimal 
approach for efficiently handling problem {2j. 

1 The exhaustive search method in (6] samples the achievable rate region 
1Z by discretizing the cross-link interference into a finite number of levels. 
Let M be the number of discretization levels. This method then needs to list 
a total number of M K ^ K ~^ rate tuples, and finds the one with maximum 
SiLi a iRi- F° r a rough case of M = 10 and K = 3, this method requires 
to search over 10 B rate tuples, which is computationally prohibitive. 



3. PROPOSED CONVEX APPROXIMATION METHOD 

3.1. Closed-Form Expression of Outage Probability 

While the probability constraints in J2bt seem intractable, there ac- 
tually exist closed-form expressions. To show this, it is noted that 
each of the probability in \2b\ can be expressed as 

luff 2 
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which is the left tail probability of the ratio of the exponential ran- 
dom variable | hfi Wi | 2 to the sum of independent exponential ran- 
dom variables |h^iUfe| 2 for k 7^ i. According to (§] Appendix I], 
l[3j has a closed-form expression as 
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Hence problem (O can be equivalently represented by 
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where pi = 1 — e,. It can be seen that ([5]l is a nonconvex optimiza- 
tion problem. Next, we show how to approximate problem (J5} by a 
convex optimization problem. 

3.2. Proposed Conservative Formulation 

The approximation method to be presented is conservative, in the 
sense that the obtained approximate solution is guaranteed to be fea- 
sible to problem (2). To illustrate the proposed method, let us define 
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where Xki,yi,Zi G R are introduced slack variables for k,i = 
1, . . . , K, and Tr(-) denotes the trace of a matrix. Substituting ([6j 
into yields the following problem 



max y^ ctiRi 
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where K\ = {1, . . . , K}\{i}, and {7h} is due to «[6c]>. Notice that 
we have replaced the equalities in doal l and j6b\ with inequalities as 
in f7cl to f7fb . It is not difficult to verify that all the inequalities in 
f7cl to l [7ft would hold with equalities at the optimum; otherwise a 
larger optimal weighted sum rate can always be obtained. Therefore, 
problem (0 is equivalent to problem l[5j. 

One can see that the objective function and most of the con- 
straints of problem ((7} are convex, except the constraints in i fTcl l 
and J7et . and the nonconvex rank-one constraints in J7ht . Let 
({■Wi}^ =1 , {Ri\f =1 ) be a feasible point of problem l[2j- Define 



x k i = \n(w"Q ki Wk), k € JCf, 



(8a) 
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for i = 1, . . . , K. Then {{xki}k^i, Vi)f=\ together with Ri, xu = 
]n(wfQiiWi), Wj = Wiwf and Zi = e Si " s " for i = 1, . . . , K, 
are feasible to problem 0. We aim to conservatively approximate 
fTct and (|7e) with respective to the point {{xki} kl a, iji}fLi- Since 
e Xki is convex, its first-order approximation at x k i, i.e., e Xki (x k i — 
Xki + 1), is a global underestimate of e Xki . Hence it is sufficient to 
achieve ( fTcl by considering the following linear constraint 



Tr(W fc Q fcl ) < e Xk '{x kl - x kl + 1), 



(9) 



for k G ICi . To approximate < f7et , we consider the following lower 
bound for e Vi + 1: 
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one. Surprisingly, it is found in our simulations that problem J 1 2b 
always yields rank-one optimal {Wj}^, i.e., Wi = Wi(wi) H for 
all i, provided that Wi 7^ 0. This implies that, for the problem 
instances we tested in simulations, an approximate beamforming so- 
lution to problem (O can be directly obtained by decomposing the 
optimal {Wj}*] of problem d!2t . 

3.3. Sequential Conservative Approximations 

The formulation l !12t is obtained by conservatively approximating 
problem ([2J with respect to the feasible point ({wi}('L 1 , {Ri\f =1 ) 
[see (8}]. It is possible to further improve the approximation per- 
formance by solving problem d!2t iteratively with the optimal 
({-Wi}fL 1 , {7?i}fL x ) at the current iteration used as the feasible 
point ({wi}fL ly {Ri\^ =1 ) for the next iteration. The proposed se- 
quential approximation algorithm is summarized in the following 
Algorithm 1: 

Algorithm 1 Proposed sequential conservative approximation algo- 
rithm for solving problem {2]l 

1: Input a feasible point ({wi}^L 1 , {Ri}^Lj) of problem (0, and 

a solution accuracy 8 > 0. 
2: Obtain {{x k i}k^,yi}?=i by {8} and obtain = 

e Vi/(eVi + i) and i2 (^) = 1 /( e s> + 1) fori = 1, . . . , if. 
3: Solve problem d 12b to obtain the optimal beamforming matrices 

{Wt}f =1 and rates {Rt}? =1 . 
4: Obtain w* by decomposition of W* = w*(w*) H for i = 

1,...,K. 

5: Output the approximate beamforming solution (wl , . . . , w* K ) 



and achievable rate tuple (R\ 



,R*k) lf I E»li a iRi 



where Oafa) = e s */(e s * + 1) and i2 (&) = l/(e Si + 1). Equa- 
tion d 1 Ob is obtained from the inequality of arithmetic and geometric 
means. By 1101) . a sufficient condition for f7et can be obtained as 



Yl,f=i a iRi\/Y^,f=i < <^ otherwise update Wi := w* 

and Ri := R* for all i, and go to Step 2. 

A feasible point to initialize Algorithm 1 can be easily obtained 
by some heuristic transmission strategies. For example, one can ob- 
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for i — 1, . . . , K, which are convex constraints. By replacing f7cl > 
and dTel l with ^ and i ll It . respectively, and by ignoring the noncon- 
vex rank-one constraints in J7ht . we obtain the following conserva- 
tive formulation for problem Q: 
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simple maximum-ratio transmission (MRT) strategy. In this strategy, 

/p<7i 



the beamforming vectors {tVi} i=1 are simply set to Wi 
where qt G C^*, ||qi|| = 1, is the principal eigenvector of Qa for 
i = 1, . . . ,K. For the ith transmitter-receiver pair, the associated 
ti -outage achievable rate of MRT is given by the maximum Ri that 
satisfies the following inequality [see $5$] 
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where 8«(&) = (^(yO)"* 1 ' 170 ^^)) 9 - 2 ^ . Problem 03 is a 
convex optimization problem; it can be efficiently solved by standard 
convex solvers such as CVX (7). 

The idea of removing the nonconvex rank-one constraints of 
{Wilf?-! in d 12b is known as semidefinite relaxation (SDR) in con- 
vex optimization theory f9]. SDR is in general an approximation 
because the optimal {Wi}|Li of problem i 1 2t may not be of rank 



Analogously, one can also obtain a feasible point of l[2j by the zero- 
forcing (ZF) transmission strategy, provided that N t is sufficiently 
large. In the next section, we present some simulation results to 
demonstrate the efficacy of the proposed approximation algorithm. 

4. SIMULATION RESULTS AND DISCUSSIONS 

In the simulations, we consider the multiuser MISO-IFC as de- 
scribed in Section 2. For simplicity, all the receivers are assumed 
to have the same noise power, i.e., a\ — ■ ■ ■ — a 2 K = a 2 , and all 
the power constraints are set to one, i.e., Pi = • ■ • = Pk = 1. 
The channel covariance matrices Qki were randomly generated. We 
normalize the maximum eigenvalue of Q«, i.e., A max (Qii), to one 
for all i, and normalize A max (Qfc,) to a value 77 G (0, 1] for all 
k G K,\, i = 1, . . . , K. The parameter 77, thereby, represents the 
relative cross-link interference level. If not mentioned specifically, 
the ranks of Qki are all set to Nt. We consider the sum rate maxi- 
mization problem by setting «! = ■■■ = or- = 1 for problem l[2}, 
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Fig. 1. Average achievable sum rate versus r\ for K = 2, JVt = 4, 
and rank(Qfci) = 4 for all fc, i. 

and set ei = • • • = eje = 0.1, i.e., 10% outage probability. For the 
proposed approximation algorithm (Algorithm 1), we set 8 = 10~ 2 
and use CVX (7) to handle the associated problem d 1 2b . All the 
simulation results were obtained by averaging over 500 trials. 

In the first example, we examine the approximation accuracy 
of the proposed method by comparing with the optimal sum rate 
obtained by the exhaustive search method in |6j. Figure 1 shows the 
simulation results of average achievable rate versus r\ for K = 2 
and N t — 4. The achievable rate of the simple TDMA scheme 
is also shown in this figure. Firstly, one can see from this figure 
that the sum rate achieved by the proposed method approaches that 
of TDMA with increased r\\ TDMA exhibits a constant sum rate 
for all r\ because there is no cross-link interference for this scheme. 
Secondly, we observe that the proposed method can exactly attain the 
average optimal sum rate for l/cr 2 = dB and l/cr 2 = 10 dB. For 
l/cr 2 = 20 dB and for r\ > 0.5 (interference dominated scenarios), 
it can be observed that there is a small gap between the rate achieved 
by the proposed method and the optimal rate. Nevertheless, this gap 
is within 3% of the optimal sum rate on average. 

In the second example, we compare the proposed method with 
the MRT scheme and TDMA for N t = K = 4. Figure 2 shows 
the results of average sum rate versus l/cr 2 . Note that, for the case 
of K = 4, the exhaustive search method in (6) is too complex to 
implement, and thus no result for the optimal sum rate is shown. 
From Fig. 2, we can observe that the proposed method achieves the 
highest sum rate among the three methods, no matter when 77 = 0.2 
or 77 = 1. One can also see that, for 77 = 0.2 and l/cr < 5 dB, 
MRT can yield a sum rate comparable to the proposed method and 
outperforms TDMA; whereas TDMA performs better for 77 = 1. 

In order to compare with the ZF scheme, in the third example, 
we extend the number of antennas to 8 (Nt = 8) and constrain the 
ranks of all channel covariance matrices to 2. The simulation results 
are shown in Fig. 3. As seen from this figure, the proposed method 
still performs best compared to the other three schemes. On the other 
hand, one can see that ZF can achieve a higher average sum rate than 
TDMA, and also outperforms MRT for 77 = 1. 
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